
iniQiuiic^iul Ar{hiv« of tilt Mil vdl Poii^roduiit School 

Calhoun: The NPS Institutional Archive 
□Space Repository 



Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


2004-09 

Analysis of the influence of turbulence and 
environmental variability on broadband 
acoustic coherence 

Eroglu, Ozer 

Monterey, California. Naval Postgraduate School 

http://hdl.handle.net/10945/1436 

Copyright is reserved by the copyright owner. 

Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


Caflwuo is the Naval Postgraduate School's public access digital repository for 
research mate rials and institutiional publicatkins created by the NPS community. 
Calhoun is named for Professor of Mathematics Guy K. Caftiouo, NPS's first 
appointed — and published — schoteily author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 Univefsity Circle 
Monterey, California USA 93943 


htt p://w ww. n ps. e du/l ib ra ry 







NAVAL 

POSTGRADUATE 

SCHOOL 

MONTEREY, CALIFORNIA 

THESIS 


ANALYSIS OF THE INFLUENCE OF TURBULENCE 

AND ENVIRONMENTAL VARIABILITY ON 

BROADBAND ACOUSTIC COHERENCE 

by 


Ozer Eroglu 


September 2004 


Thesis Advisor; 

Kevin B. Smith 

Thesis Co-Advisor; 

Samuel E. Buttrey 


Approved for public release; distribution is unlimited. 




THIS PAGE INTENTIONALLY LEET BLANK 



REPORT DOCUMENTATION PAGE 


Form Approved 0MB No. 
0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including 
the time for reviewing instruction, searching existing data sources, gathering and maintaining the data needed, and 
completing and reviewing the collection of information. Send comments regarding this burden estimate or any 
other aspect of this collection of information, including suggestions for reducing this burden, to Washington 
headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 
1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project 
(0704-0188) Washington DC 20503. 

1. AGENCY USE ONLY (Leave 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

blank) September 2004 Master’s Thesis 

4. TITLE AND SUBTITLE: Title (Mix case letters) 

Analysis Of The Influence Of Turbulence And Environmental Variability On 

Iroadband Acoustic Coherence _ 

6. AUTHOR(S) LtJG Ozer Eroglu_ 


11. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect the official 
policy or position of the Department of Defense or the U.S. Government. 


13. ABSTRACT (maximum 200 words) 

Acoustic propagation in the littoral regions of the world, even over short ranges, can be complex at high 
frequencies (> 5 kHz), and applications such as underwater detection and communications suffer as a result. To 
this end, the Asian Seas International Acoustics Experiment (ASIAEX) was conducted with funding from the 
Office of Naval Research. One phase of this experiment took place from 29 May to 9 June 2001 and focused on 
short-range, shallow water acoustic propagation in the East China Sea. This thesis will be based on some of the 
measurements taken during the East China Sea experiment by a team from the Applied Physics Laboratory of the 
University of Washington. The environmental parameters and array geometry used in this and previous modeling 
work are as close as possible to the experiment. 

The objective of this thesis is to better understand the nature of turbulent perturbations and how they affect 
short-range acoustic propagation in a shallow water environment. This will include variations in total energy of 
the turbulence, variations in associated length scales, and variations in depth. In addition, the influence of other 
variable factors on signal coherence will also be examined. Most notably, variations in background sound-speed 
structure and source depth variability will be included._ 


16. PRICE CODE 


NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. 239-18 


20. LIMITATION 
OF ABSTRACT 


15. NUMBER OF 
PAGES 95 


14. SUBJECT TERMS 

Shallow water variability, turbulence, vertical coherence, sound speed curvature 


18. SECURITY 
CLASSIFICATION OF THIS 
PAGE 

Unclassified 


19. SECURITY 
CLASSIFICATION OF 
ABSTRACT 

Unclassified 


17. SECURITY 
CLASSIFICATION OF 
REPORT 

Unclassified 


12b. DISTRIBUTION CODE 


12a. DISTRIBUTION / AVAILABILITY STATEMENT 

Approved for public release; distribution is unlimited. 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Naval Postgraduate School 

9. SPONSORING /MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

N/A 


5. FUNDING NUMBERS 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


1 



























THIS PAGE INTENTIONALLY LEET BLANK 


11 



Approved for public release; distribution is unlimited. 

ANALYSIS OF THE INFLUENCE OF TURBULENCE AND ENVIRONMENTAL 
VARIABILITY ON BROADBAND ACOUSTIC COHERENCE 

Ozer Eroglu 

Lieutenant Junior Grade, Turkish Navy 
B.S., Turkish Naval Aeademy, 1998 

Submitted in partial fulfillment of the 
requirements for the degrees of 

MASTER OF SCIENCE IN ENGINEERING ACOUSTICS 

& 

MASTER OF SCIENCE IN OPERATIONS RESEARCH 

from the 


NAVAL POSTGRADUATE SCHOOL 
September 2004 


Author: Ozer Eroglu 

Approved by: Kevin B. Smith 

Thesis Advisor 


Samuel E. Buttrey 
Co-Advisor 


Kevin B. Smith, Chairman 

Engineering Aeousties Aeademie Committee 

James N. Eagle, Chairman 
Department of Operations Researeh 

iii 



THIS PAGE INTENTIONALLY LEET BLANK 


IV 



ABSTRACT 


Acoustic propagation in the littoral regions of the world, even over short ranges, 
ean be eomplex at high frequeneies (> 5 kHz), and applieations sueh as underwater 
deteetion and eommunieations suffer as a result. To this end, the Asian Seas International 
Aeousties Experiment (ASIAEX) was condueted with funding from the Offiee of Naval 
Researeh. One phase of this experiment took plaee from 29 May to 9 June 2001 and 
foeused on short-range, shallow water aeoustie propagation in the East China Sea. This 
thesis will be based on some of the measurements taken during the East China Sea 
experiment by a team from the Applied Physies Eaboratory of the University of 
Washington. The environmental parameters and array geometry used in this and previous 
modeling work are as close as possible to the experiment. 

The objeetive of this thesis is to better understand the nature of turbulent 
perturbations and how they affeet short-range aeoustie propagation in a shallow water 
environment. This will inelude variations in total energy of the turbulenee, variations in 
assoeiated length seales, and variations in depth. In addition, the influence of other 
variable factors on signal coherence will also be examined. Most notably, variations in 
background sound-speed structure and souree depth variability will be ineluded. 
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1. INTRODUCTION 

Having the ability to perform operations in the littoral regions of the world is one 
of the most erueial issues if a navy intends to be fully eapable. Sinee these regions differ 
from the so-ealled “blue waters” in many ways, even the most thorough deep water 
research may not be applicable. One of the most extensive researches conducted in 
littoral regions was the Asian Seas International Acoustics Experiment (ASIAEX), 
funded by the Office of Naval Research. The second phase of this project took place in 
the East China Sea (ECS). 

One component of the ECS experiment involved the analysis of acoustic signal 
coherence over the length of a vertical array at relatively short range (500 m) from a 
controlled source. Eor certain source/receiver geometries, direct path arrivals, single 
surface bounce, and single bottom bounce paths could be distinguished, allowing analysis 
of coherence for those paths independently. 

Previous thesis work done by ET Stephen C. Karpi at the Naval Postgraduate 
School contributed to the analysis of the direct path (no boundary reflections) arrivals by 
numerically evaluating the influence of various shallow water sound speed perturbations. 
Karpi’s research employed vertical coherence to measure the effects of longer scale, 
internal wave-like structures and smaller scale, turbulent-like fluctuations on the direct 
path propagation. His work showed that perturbations by an internal wave-like structure 
affected the coherence only slightly. However, turbulent-like fluctuations caused 
significant signal decorrelation RMS. 

This thesis examines both the model geometry, including source and receiver 
depths, and the environmental variables, which are background sound-speed profile, 
turbulence strength and the turbulence outer length scale. The general geometry and 
source characteristics used during the East China Sea experiment by Peter Dahl of the 
Applied Physics Eaboratory of the University of Washington are the basis of numerical 
studies in this thesis. 
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The model used in this work is the Monterey-Miami Parabolie Equation (MMPE) 
propagation model. The MMPE model was developed by Smith and Tappert in 1996. The 
speeifie version used in this work ineorporates the range-dependent turbulent fluetuations 
and the souree speetrum eharaeteristies employed in the ASIAEX experiment. This 
model provides an effieient algorithm for eomputing the broadband souree response in 
sueh range-dependent environments. 

A realistic turbulence spectrum has multiple parameters which depend on various 
features of the ocean environment. In this work, only two parameters, a strength 
parameter to adjust the overall RMS value, and an outer length scale beyond which 
fluctuations have no contribution, will be considered variables. As a reference, 10 m is 
used for the outer length scale, and four other different values are tested. Because only 
weak scattering is expected, the outer length scale is not expected to directly affect the 
propagation. However, the dominant contribution to the scattering is expected to be at the 
Eresnel radius, and this is on the same order as the reference scale of 10 m. Thus, we may 
expect to see some differences with changes in the outer scale due to the energy in the 
turbulent perturbation near the Eresnel scale. Conversely, the turbulence strength 
parameter is expected to have a direct effect on scattering, and thus affect the signal 
coherence. 

Besides the two aforementioned parameters of turbulence, the effects of variable 
background sound-speed profiles are examined in this research. The reference sound- 
speed profile (SSP) is taken as the average value of the measurements from ASIAEX. 
Three more SSPs are also examined. One of them is very close to the reference case, but 
has more variability along the depth. The other two are very different from the reference 
case. One of them has a much deeper thermocline and steeper gradient below, while the 
other has a shallower thermocline and generally less steep gradient below. The effects of 
SSP are examined using not only the vertical correlation but also using some simple ray 
tracing analysis involving sound-speed relative curvature (Duda et al, 1994). 

The last variables investigated in this work are source and receiver depths. Since 
the model computes the propagation at all depths, no additional computational effort is 
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required to obtain the solution at any reeeiver depth. To eompute vertieal eoherenee, 
short sub-arrays are modeled by extraeting the solution over 5 m vertieal spans. The 
reeeiver sub-array depths that are analyzed are 35-40, 45-50, 55-60, and 65-70 m. The 
referenee value used for souree depth is 50 m, while 30, 40, 60, and 70 m are also tested. 
It should be noted that above 30 m or below 70 m of the souree or reeeiver depths, the 
direet path propagation will overlap with those arrivals from surfaee seattering and 
bottom refleetions, for the 2 msee pulse modeled. Sinee the foeus of this work is to deal 
with the direet path propagation, these depths are disregarded in the analysis range. 

Chapter II eontains an overview of the theory behind the MMPE model and the 
turbulent speetrum. The implementation of these models and the post-proeessing methods 
are also examined in this ehapter. Different model parameters and their range of analysis 
are eovered in Chapter III. Chapter IV foeuses on post-proeessing and the analysis of 
modeled data. The last ehapter presents eonelusions derived from the analysis and 
suggestions for future researeh. 
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II. NUMERICAL METHODS AND IMPLEMENTATION 


A. MONTEREY-MIAMI PARABOLIC EQUATION (MMPE) MODEL 


The parabolic equation (PE) method was introduced into underwater acoustics in 
the early 1970’s by Tappert. Due to its efficient numerical algorithm for computing 
solutions to general range-dependent problems, it has become the most popular method 
for computing the acoustic field in complex underwater acoustic environments. The 
MMPE model, developed by Smith and Tappert in the early 1990’s, is based upon the 
parabolic equation method. In this chapter, a simple derivation of the MMPE model is 
described by borrowing heavily from Smith’s descriptions. 

We begin by defining the time-harmonic acoustic field in a cylindrical coordinate 
system by 

P{r,z,(p,(ot) = p{r,z,(p)e '“’‘ . (1) 


Substituting Eq. (1) into the linearized acoustic wave equation in cylindrical coordinates 
leads to the Helmholtz equation. 


lA 

r dr 


dp\ 1 d^p d^p 


dr 


+ 


r d(p dz 


+ ^^^ + kln^{r,z,(p)p = -A7i;PQd{x-Xs), 


( 2 ) 


where 


= colC q , 

is the reference wavenumber, and 


( 3 ) 


n{r,z,(p) = (4) 

c{r,z,(p) 

is the acoustic index of refraction. In Eq. (4), Cq is the reference sound-speed and 
c{r,z,(p) is the acoustic sound-speed. Note that c{r,z,(p) carries all the information 
about the environment, except the density, which is unimportant in this development 
since we are only interested in propagation in the water column. Eq. (2) implies that the 
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The last term in Eq. (6) deereases quickly due to 1/r^ dependence and is generally 
neglected. The second term accounts for azimuthal coupling between different radials 
and is small enough to consider as negligible. Neglecting this term is referred to as the 
“uncoupled azimuth (UNCA) approximation.” 

The remaining Helmholtz equation can be further simplified by introducing the 
operator notations 

and Q^^={^ + s + \y , (7) 

where 


2 1 15 ^ 


Eq. (6) can then be represented in the form 


(F +Fap>=o. 

Eactoring Eq. (9), outward and inward fields are represented in the form, 
(Pop + iKQop)(Pop - iKQop > + iK{Pop . Qop ]w = 0 


( 8 ) 

( 9 ) 

( 10 ) 


Assuming weak range dependence, the commutator \_Pgp,Qgp\ is negligible. 
Eurthermore, defining the outward propagating field as 



(11) 
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the outgoing solution satisfies 

PVT/ 

Po,^ = iKQo,^ or = (12) 

If the baekseattered energy is negligible, Eq. (12) represents the eomplete description of 
the forward propagating acoustic energy in the waveguide. This equation is the 
foundation for all underwater acoustic PE models. 

To solve the parabolic equation given in Eq. (12), first the field should be 
decomposed into two parts, an envelope function and a phase term. A reasonable 
decomposition for the PE field function is defined by 

‘P = . (13) 

Using this modulating function, we get the PE for the field function in the form 

■ ( 14 ) 

PE solutions are based on developing an approximation to the pseudo-differential 
operator . There are basically three methods to solve the PE: the Einite Element (EE) 

method; the Einite Difference (ED) method; and the Split-Step Eourier (SSE) method. An 
approximation well suited for the SSE algorithm, on which the MMPE model is based, 
was introduced by Thompson and Chapman (1983) and is commonly referred to as the 
wide-angle PE (T-C WAPE). This approximation generally provides good accuracy in 
phase. 

To obtain the SSE algorithm, let us rewrite the PE in terms of the Hamiltonian- 
like operator as 

^ = -iKH.,¥ , ( 15 ) 

dr 

where 

(16) 
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This Hamiltonian-like operator defines the evolution of the PE field funetion in range. 
Solutions of the 'F funetion at different ranges ean be defined by a marehing algorithm 
of the form 


^(r + Ar) = 0(r)^(r) , (17) 

where 0(r) is a unitary operator that progresses the solution out in range. The propagator 
0(r) ean be approximated as 

, (18) 

where 


r-fAr 

HoAr) = — J . (19) 

r 

The advantage of the SSF algorithm is based on the faet that is separated into 
a combination of spatial and differential operators, which can be expressed as 

( 20 ) 

where 


and 


Uo,=- 


\l\ + s - \ 


= -{\-n) 


( 21 ) 




1 + 


kl dz^ 


( 22 ) 


Keeping in mind the assumptions made thus far, the propagator can be made 2 
order accurate in range step (Ar) by employing the form. 


nd 


0(r) = 


Ar 



Ar 

e ^ 


-i 

e 


_ 


_ 

_ 


(23) 


The operator Uop is a multiplication operator in the z-space and, hence, is a diagonal 
matrix. The operator Top is not diagonal in the z-space, but it is in wavenumber space; 

therefore, the kz-domain operator ^op is desired. 
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It is now possible to outline how the PE/SSF algorithm steps out in range. To 
begin, the PE field funetion is defined at some range, r, in the z-domain. This funetion is 
then multiplied by the operator, e defined at the beginning of the range step. 

A transformation to the ^^-domain then oceurs so that we can multiply by e . This 
result is then transferred back to the z-domain to be multiplied by the z-space operator 
defined at the end of the range step, e discrete fast Fourier Transform 

(EFT) is being used in the code assuming the convention 


y,{z) = FFT{ifr{k^)) 
and 


(24) 


y/{K) = IFFT{y,{z)) . 
The PE/SSF implementation can then be represented by 


(25) 


Ar,z) 


= e 


Af 

-ik^ — (r+Ar,z) 

”2 VfT 



-ik„^U„p(r,z) 

e ^ y/\ 


g i-KAFiip (k ^) ^ Jp'p'J' 

[r,z) 


V 

) _ 


(26) 


where 





Uoj 


(27) 


B. TURBULENT SOUND-SPEED PERTURBATION THEORY 

Acoustic propagation through a medium is subject to the properties of that 
medium. When a spherical wavefront from a point source traverses a medium containing 
fluctuations in sound speed, portions of the front travel faster or slower than other parts 
of the wavefront and the front distorts as a result. In the ocean, such sound-speed 
inhomogeneities result from phenomena like layering, internal waves, and turbulence 
(Duda, et al, 1988). In this work, it is assumed that perturbations of sound speed are 
generated purely by turbulent fluctuations. These fluctuations are characterized using 
statistical relationships developed in previous work (Tatarskii, 1961; Duda, et al, 1988). 
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The statistical theory of turbulence characterizes the turbulent perturbation field as 
realizations of a random function. Such a characterization is only valid when changes in 
signal amplitude, due to fluctuations in the medium, do not exceed the average amplitude 
of the signal. This is referred to as the “weak scattering regime” (Duda et al, 1988). 

If we define the index of refraction perturbation as 

= = ^{r), (28) 

c{r) 

then the variance of /u{r) is related to the spectral density of the perturbations, S^ {K ), 
by 


var(//(F)) = 


1 


\s^{K)d^K, 


(29) 


where the integral is a three-dimensional integral over all wavenumber space. The 
variance is defined by 


so 


var(//(F)) = (|//(F)n = —- \ //{r)d^r , (30) 

' ' volume 




volume 


\s^{K)d^K. 


If we define the spectrum (transform) of ju{r ) as FJ^JC) by 


//(r) = 


f F(K)/^*^d^K, 

{Inf I 


then the correlation of jLi{r) satisfies 


(31) 


(32) 


00 1 °° _ 2 - 

J MF)/(F-FVV==— y J d^K\F^{K)\ . (33) 

_00 ^ ’ —00 

I “ 1 ^ 

Thus, \F^{K)\ is the Fourier transform of the correlation function. It is easy to see that 
by setting r' = 0 , the square of the signal spectrum and the spectral density are related by 


10 



It follows then that 


= {volume)*S^{K) . 


( 34 ) 


\F^{K)\ = \^olume* S^{K)~^ . (35) 

Note that //(F)is unitless while F^{K) and (^)both have units of m^. 

Similar reasoning holds for a two-dimensional shoe at a fixed value of one 
variable, e.g. /ii{x,y = 0,z), sueh that 

-« < 00 00 

- - JJ 7 = 0, z)dxdz = J J {K^ , )dKJK^, (36) 

where {K^■>F^) = V^{k^,k^,ky =0) is the two-dimensional speetral density along the 

slieey = 0. Furthermore, 

00 00 

Mix, z) = —^ J J i^x, , (37) 

V '^ ^ 7 —00 —00 

where G^{k^,k^)is the speetrum (transform) of the two-dimensional perturbation field. It 
follows that 

\G/k,.k,)\^[Area»r^{k,.k,)]"\ (38) 

The two-dimensional speetral density ean be obtained from the three-dimensional 
speetral density aeeording to (Tatarskii, 1961), 

1 

rjKX)=—ls^(k,,k,,k,)dk,. (39) 

-00 

At this point in the turbulent analysis, perturbations in sound speed have been 
related to the speetral density of the turbulent fiuetuations. Developing an expression for 
the turbulent perturbation field and ealeulating its statisties remain to be done. Tatarskii 
approximates the turbulent field as isotropie and homogeneous by using the method of 
strueture funetions (Tatarskii, 1961). Aeoepting that large-seale inhomegeneities may be 
present, strueture funetions restore homogeneity by redueing the seale of the analysis. 
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Standard statistical relations are then applied, resulting in a three-dimensional spectral 
density function of the form (Tatarskii, 1961), 

E{K) = A ^ , (40) 

where , A is a sealar multiplier, and e is the energy dissipation rate. 


In the literature (e.g., Duda, et al, 1988), this speetral density has been adapted as 
follows. To avoid the speetral density beeoming infinite at if = 0, a wavenumber 
threshold (outer scale) limit is defined as kt such that 


EiK)^ 


1 


(K^+k^) 


2\ll/6 


(41) 


thus maintaining the K power law from Tatarskii’s speetrum. A highpass filter of 
the form 


is also added to force the spectral density to zero at K=0. 


(42) 


For high wavenumbers, a high-frequeney eut-off related to the Batchelor 
spectrum is given by 


where 


Rb(K) = qxp 




(43) 


A ^ A 


1/4 


kg — 


j 


(44) 


is the Batchelor wavenumber, e is the depth averaged kinetic energy dissipation rate in 
W/kg, V = 1.40e-6 m /s is the kinematie viseosity, and Kj ~1.00e-7 m /s is the thermal 
diffusivity of sea water. The eonstant q is an order unity faetor, and 3.7 is used here in 
concordance with earlier work (Duda, et al, 1988). 
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The depth averaged kinetic energy dissipation is defined as (Henyey et al, 1986), 


e=0.3 —( 



2 

cosh ' 

r N{z))\ 

\L^ ' 


K fl jj 


(45) 


The brackets ( indicate a depth average, 7 , = 1 is the characteristic mode for shallow 
water, bE^j^ = 0.5m is a measure of internal wave intensity, fj is the inertial frequency 
in rad/sec, and the buoyancy frequency N{z) may be computed from 


N{z) = g^'^ 


1 




dz 


1/2 


(rad/ sec). 


(46) 


where p{z) and (z) are the density and potential density, respectively, as a function 

of depth. This expression for "e is valid for turbulence created by internal wave decay 
(Henyey et al, 1986). It is an expression developed for deep water and ignores other 
possible sources of turbulence. 


By combining these quantities, the three-dimensional spectral density has the 

form 


(b 

S(K) = — , -v/?«(^) • 

{K^+kf "" 


(47) 


From this expression, the two-dimensional spectral density is calculated below. The 
integration is applied only to the portion defining the turbulent spectrum (the scaling 
factor and the wavenumber roll-off terms are added afterwards). 


®3Z2(^) 


(48) 


The two-dimensional spectral density is then 


j CO 

^2DiK^kJ = — \ (E,^{k^,k^,k^)dky 


1 

J-f 


I 


27riik:+k +k:+k^y 


rdk,. 
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( 49 ) 


1 V^r(i) 

27r[kl+kl+k^Y" m) 

Employing the same highpass filter, high-wavenumber Batehelor speetrum eutoff, 
and turbulent strength parameter sealing, the two-dimensional speetral density of 
turbulent fluetuations of the index of refraetion is 


VAK) = 


r(l) 4 




In r(^) K^+kf 


R,(K). 


Note that V (K) has units of m^ as expeeted. 


(50) 


C. IMPLEMENTATION METHODS 

I. Numerical Implementation of Turbulence 

In the above diseussion, V^(K), the two-dimensional speetral density of 
refraetive index perturbations, was related to /j(x, z ), the perturbation field. A random 
realization of the two-dimensional turbulent field of sound-speed perturbations is 
obtained by employing the two-dimensional speetrum, {k^ ,k^) , and defining 

.. 00 00 

(2;z-) J 

where B{k^,k^) and 6{k^,k^) are two-dimensional random number funetions. Sinee 

the eomplex amplitude of eaeh wavenumber eomponent, Be^, should exhibit a normal 
distribution, the random phase of eaeh eomponent is a uniformly distributed random 
variable on the unit eirele. In other words, we define 

6 = 2nr ^, (52) 

where is a uniformly distributed random variable in the interval [0,1]. The magnitude 
B, however, exhibits a Rayleigh distribution. Because we are in fact generating a 
realization of the spectral density, we must consider the magnitude-squared, which has a 
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negative exponential distribution. Thus, the random amplitude of each wavenumber 
component is obtained by 


B = (53) 

where Tj is another independent, uniformly distributed random variable in the interval 

[ 0 , 1 ]. 


Since we are implementing this in a model with a discretely sampled grid of the 
environment, we must consider discrete Fourier transforms computed using an FFT 
algorithm. Thus, with scaling appropriate to the FFT algorithm being employed. 


M{X^,Z,) = 




M N 

Z Z [g, (*„. K. ■ Kn ’ ] 

m=\ n=\ 


e 


(54) 


From the previous analysis of the relationship between the perturbation spectrum and the 
spectral density. 


(55) 

where we here define the spectral density up to a scaling factor 


UK) = 


r(l) i 




R,(K). 


-2 


(56) 


Itt r(^) (K^+k^y'^ K^+kf 

Note that is the fixed turbulent scaling parameter and may be adjusted to produce the 
desired RMS perturbation. 

Combining the above, we find 


AjL _ M N 


i(K, 


-KnH) 


{In) 


n=\ n=\ 


= y/AreaFFT^° (K)B(K)e‘^^^^ ] (57) 

For a given computational area, the factor Area = RZ, where the turbulent field exists. 
Furthermore, 
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( 58 ) 


so 


Thus 


^k^^ky 


Area 


Area = 


Ak^Ak 

A y 


(59) 


M{X.,Z,) = ^Ai^FFT^^ [rj\K)B{K)e^‘^^‘^^\ (60) 

Once a realization of //(x, z) is computed, we simply multiply by Co to obtain the 
sound-speed perturbation, 

5c{x,z) = c^iu{x,z). (61) 

Teehnioally, we might multiply by c(x,z), a sound-speed profile with range. However, 
variations in c(x,z) are typically less than two percent and so using Co is a reasonable 
approximation. Furthermore, note that we are employing the perturbation in eylindrieal 
eoordinates, replaeing x by r, so 

Sc{x,z) ^ Sc{r,z). (62) 

This is also justified so long as we only eonsider propagation along the two-dimensional 

shoe in (r,z). See Figure 1 for an example realization of the sound-speed perturbation 
field. 
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Realization of Turbulent Sound Speed Perturbation 
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Figure 1. Sample realization of the turbulent sound speed perturbations in m/s. The RMS 
sound speed fluctuation is 0.41 m/s. 

To compute the parameters A^(z) and s , we need to extract both the buoyancy 
frequency, N(z), and the variance of the temperature gradient, ^ The former 

may be defined as 

1 dp(zy'^" 


N{z) = g 


1/2 


(rad/ sec), 


(63) 


p{z) dz 

where p{z) and Pp{z) are the depth profiles of the density and potential density, 

respectively, and g is the acceleration of gravity. Both densities were measured during 
the ASIAEX experiment by CTD casts. 
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2 . 


Creation of Source Spectrum 


During the experimental periods of interest, characteristic continuous wave (CW) 
and frequency modulated (FM) pulses were transmitted. The CW pulses of 2 to 3 msec in 
duration had center frequencies of 4, 8, 16, and 20 kHz. The FM pulses are not analyzed 
here. To simulate a source pulse of 2 msec, an analysis must be performed in the time 
domain after the model has run. A band-pass filter of approximately 8 kHz must also be 
incorporated, to be consistent with the work being done at APL-UW. However, before 
these concerns can be addressed in detail, a brief summary of time domain processing 
within the MMPE is necessary. 

Recall from the development of the parabolic equation that the time-harmonic 
acoustic field (ignoring azimuthal dependence) is given by 

P{r,z,o)t) = (64) 

A discrete FFT algorithm is used to represent the field in the time domain as 

P{r,z,t) = FFT[p^{r,z)\ = Y,pSr,z)e '‘“‘ ■ (65) 

CO 

In other words, Eq. (64) represents a single chosen frequency and a single component of 
the time dependent field. Thus, for analysis in the time domain, the model must be run at 
multiple frequencies. The inputs required by the model for such analyses are a center 
frequency ifo), a bandwidth {BW), and the number of frequency bins to compute 
(transform size, N). In this research, the bandwidth used was 8184 Hz with a transform 
size of 1024, consistent with Karpi’s previous work. 

To view the arrival structure at the receive array (r = R, whereR is the range of 
the array from the source), we must compute the complex field p^^iz) for many 

frequencies and then Eourier transform to obtain Pj^{z,t) , the complex pressure values in 
time/depth space. An example P,^{z,t) is shown in Eigure 2. Outside the given bandwidth 
of 
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( 66 ) 


f , 


( _ BW\ 

fc-^ 

</< 

( 

+ 

V 2 ; 


^ 2 ; 


the source amplitude is zero. The MMPE model ignores the zeroes outside the bandwidth 
and computes a base-banded field to avoid unnecessary computations. The resulting 
expression is 


where =l7tf^. The complex pressure values ^(z) are thus placed in frequency 
bins symmetrically about co-co^ =0 with limits of ±^BJV . 



Transmission Loss (dB re 1m) 
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Figure 2. Pulse arrival structure for a source depth of 50m, RMS sound speed 
perturbation of about 1 m/s, and an outer cutoff length of 10 m. 

Note the transform pair actually used is 
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( 68 ) 




icoR/ 


where is the PE field function and the overall phase factor e is 

neglected. As a result, arrival times are given as values of “reduced time”, r, such that 


r = 


^ r'^ 

t - 

V 


(69) 


where is the reference sound speed. The complex pressure values P'i^{z,t) are thus 


placed in time bins symmetrically about r = 0, with limits of ’ where T is the total 
length of the time window. The time step is given by 


At =-« .122msec , 

BW 


(70) 


with a corresponding frequency step of 


A/ = 


BW 


N-\ 


^ = 8Hz. 
1024-1 


(71) 


Running the model at multiple frequencies has thus enabled the analysis of the time 
dependent field. 


The APT team from the University of Washington attempted to remove some of 
the out-of-band noise from the ASIAEX data by applying a digital band-pass filter to the 
signal. During the filtering process, they preserved 8 kHz of bandwidth, using roughly a 
+/- 4 kHz band around the center frequency. However, the band-pass filter was not 
always applied symmetrically about the center frequency of the signal. In order for the 
model’s results to be comparable, the model was run using a center frequency defined by 
the center of the filter. For example, on a model run meant to simulate a 2 msec pulse at 
20 kHz, a center frequency of 19 kHz was input to the model, corresponding to the center 
of the filter. Thus, the output data is centered properly, but still representative of a 
continuous source. 
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To simulate a 2 msec pulse, a filter function was applied to the data. A square 
filter function (shown below in Figure 3) was chosen to represent the pulse in the time 
domain. The length of the filter corresponds to the number of discrete time steps that are 
needed to comprise the pulse, found by 


1 


0.8 


■o 0.6 

3 

Q. 

E 

< 

0.4 


0.2 


0 

-0.06 -0.04 -0.02 0 0.02 0.04 0.06 

Reduced Time (sec) 


Square filter representing 2 msec pulse 


- 2 msec 


Figure 3. Square filter function representing a 2 msec pulse. 


2 msec , , 
=-= 16, 


pulse 


At 


(72) 


where the time step is 0.122 msec, as shown in Eq. (70). This square filter will then be 
transformed to a source spectrum in the frequency domain, as shown in Figure 4, and 
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applied to the output data. However, the output data is centered at the middle of the 
bandpass fdter, and the source spectrum should be applied at the center frequency of the 
source. Thus, a shift may be required. In the above example, a shift of 1 kHz would be 
needed to center the source spectrum at 20 kHz instead of 19 kHz. Table 1 provides a list 
of the fdters used and the corresponding shift required. 


Experimental 
Carrier Freq. (kHz) 

Minimum Freq. of 
Filter (kHz) 

Maximum Freq. of 
Filter (kHz) 

Center Freq. of 
Filter (kHz) 

Required Shift 
(kHz) 

4 

2 

10 

6 

2 

8 

5.5 

13.5 

9.5 

1.5 

16 

12 

20 

16 

0 

20 

15 

23 

19 

1 


Table 1. Characteristics of the bandpass filters used in data processing for each 
experimental carrier frequency. The required shift is the difference between the center of 
the filter and the experimental carrier frequency. 


To allow for shifting, the generated source spectrum must have a wider bandwidth 
(BJV) than the output data. The bandwidth is simply doubled, resulting in a new time step. 
At ', given by 

At'= —-— = —-—« 0.061msec, (73) 

BJV' 2BW 


and thus twice as many steps are required for a 2 msec pulse. 


, 2msec 

n,uise - - 32 . 


(74) 


Note that the total time window T and the frequency step Af remain constant. The end 
result is an identical spectrum with more frequency information added on each side to 
make the spectrum twice as long. The original spectrum can be seen below in Figure 4 
and Figure 6 illustrates the extended version. 
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Source spectrum 
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Figure 4. Source spectrum created from 2 msec square window with a bandwidth of 8184 
Hz. 


In the frequency domain, the square window takes the expected sine function 
form, given by 


sinc(x) = 


sin(x) 


X 


(75) 


The source spectrum developed above must now be combined with an appropriate band¬ 
pass filter. A rectangular window with a Hanning taper at each end was chosen for this 
purpose. The Hanning window that was used can be defined as 




f 

6 k 

1 -cos 

2 ;r 

V 



where k= l,...,n. 


(76) 
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The shape of a Hanning window can be seen in Figure 5. The Hanning window 
was subsequently divided into halves, each half being applied to an end of the rectangular 
window, and the total filter width being equal to 8 kHz in this case (spanning 15-23 kHz). 
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Figure 5. Typical Hanning Signal-Processing Window 

This system response filter was then used as a multiplier on the 2 msec pulse 
spectrum, thus completing the source amplitude function. An example bandpass filter 
and source spectrum are shown together in Figure 6 for a carrier frequency of 20 kHz. 


Hanning Window 
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Source Amplitude Function 



Figure 6. Extended source spectrum overlaid with 8 kHz bandpass filter. In this case, a 
carrier frequency of 20 kHz is used, thus requiring a 1 kHz shift. 

3. Direct Path Extraction 

At the boundaries of the medium, the acoustic signal interacts with the surface of 
the water and the iso-speed bottom (set at 1600m/s). In order to remove the interface 
effects from the analysis and focus strictly on the volume effects, the direct path arrival 
was extracted. As can be seen in Figure 7, for a source depth of 50 meters, several 
propagation paths are evident in the pressure time series (received at 500m range). 
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Arrivals at Shallow Array, Lo = 10 meters, rms = 1.125m/s 


xIO’^ 




Figure 7. Pressure time series reeeived at shallow and deep depths at a range of 500 
meters. Lo is the outer eutoff length, and RMS is the RMS value of the sound speed 
perturbation. 

A source depth of 25 meters in some cases results in a direct path signal 
indistinguishable from the surface bounce arrival. As in earlier research by Karpi, a 50- 
meter source depth will be considered as a reference. Plots similar to Figure 7 were 
analyzed for each model run to determine the times encompassing the direct path arrival. 
The time series was then extracted directly, without filtering, as shown in Figure 8. The 
direct path extraction is performed at six varying depths across each of the roughly 5 m 
length vertical array clusters. The uppermost depths of each sub-array were defined as 
the numerical mesh points closest to 35, 45, 55, and 65 m. The depths of the five lower 
elements were then approximately 15, 45, 105, 225, 465 cm below the uppermost 
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element, thus spanning roughly 5 m. Therefore, 24 total signals are extracted, one for 
each of the six elements in the four sub-arrays. 



Figure 8. Extraction of the Direct Path Arrival. 

4. Vertical Spatial Correlation 

It is now possible to obtain a quantitative estimate of signal degradation by cross- 
correlating the extracted direct path signals. Correlation is a measure of the similarity 
between two waveforms. In this case, signals recorded at different depths along the same 
propagating wavefront are compared. If signal degradation was not present, such as when 
propagating through a homogeneous medium, we would expect the received pressure 
signal to be very similar at different points along the front. As the similarity of the 
received signals along the front decreases due to inhomogeneities, we say that our 
wavefront has degraded. 

In the experiment, each of the six elements in a vertical array cluster has a 
corresponding received pressure signal recorded at the depth of the element. Thus, 15 
different cross-correlations are possible (element 1 with element 2, element 2 with 
element 3 and so forth), in addition to the auto-correlation, each with a depth separation 
given by the separation of the elements being correlated. The final output is a graph of 
the normalized maximum vertical cross-correlation as a function of the depth difference 
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of the elements being correlated with 16 points. These 16 points include each possible 
combination between elements and the autocorrelation of the signal at element one. 


Practically, each extracted direct path signal is cross-correlated using the 
MATLAB function xcorr, which generates one number for the normalized max 
correlation. Note that “normalized” means the autocorrelation (the correlation of the 
signal with itself) is assigned a value of 1. The normalized cross-correlation can be 
represented as follows: 


C,,,2ir) = e 





\pi \T)p^\T -T)dT 


' A \T) 

Ut\ 

A \T) 

"dT 


(77) 


where and p^* are the first base-banded pressure signal and the complex conjugate of 
the second base-banded pressure signal, respectively. 
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III. MODEL VARIABLES 


A. EXPERIMENTAL AND MODEL GEOMETRY 

The input parameters to the MMPE model were taken mostly from data sets of the 
East China Sea experiment, and the geometry and source characteristics were modeled 
after the short-range experiment conducted by the group from APE-UW. The horizontal 
range between the source, placed at a depth of either 25 or 50 m, and receivers was set at 
500 m, and the water column depth was set at 100 m. In the ASIAEX experiment, two 
four-element vertical arrays were placed at nominal depths of 26 and 52 m. The 
distances between the elements of each array were 13, 30, and 60 cm. Eigure 9 below 
shows the experimental configuration. 
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Eigure 9. Experimental Geometry Ifom ASIAEX 
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The MMPE model is capable of incorporating both shallow and deep bottom 
characteristics. In this work, however, the deep bottom depth was set at 3000 m, deep 
enough not to affect propagation. The sound speed of the shallow bottom was set at a 
constant of 1600 m/s. This version of the MMPE model requires a reference sound speed, 
Co, set at 1500 m/s for this research. It should also be noted that all results presented here 
are for runs done in “efficiency” mode as opposed to “accuracy” mode. The difference 
between the two modes is in mesh size, Ar, Az . Efficiency mode was chosen because it 
is much faster than accuracy mode, and in most cases should be accurate enough for our 
direct path analysis. 

Though this research is mostly concerned with the effects of environmental 
variables, it is also a point of interest to learn about the effects of the geometric variables, 
such as source and receiver depths. Since the model computes the solution at all 
computational depth grid points, it is a simple matter to analyze the dependence of signal 
coherence on receiver depth. Thus, the numerical work presented here differs in 
geometry from ASIAEX in that a total number of four receiver sub-arrays with six 
elements each are processed. All analysis in this work used uppermost sub-array element 
depths of 35, 45, 55, and 65 m, with total apertures of roughly 5 m. 

1. Source Depth 

In contrast to the ASIAEX geometry, the influence of live different source depths 
is examined in this research. A depth of 50 m is assumed as the reference source depth. 
Other source depths included 30, 40, 60, and 70 m. The reference geometry for the 50 m 
source depth and the four sub-arrays is displayed in Eigure 10. 
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Figure 10. Model Geometry Used in This Research 


B. ENVIRONMENTAL MODELS 

The focus of this research is on the influence of turbulence, so parameters 
associated with the description of the turbulent field will be varied to determine their 
influence. Specifically, we will examine the effect of variable outer length cut-off scale, 
kt, and overall turbulent strength scaling parameter, These parameters have a direct 
effect on the structure of the turbulence, and thus a direct effect on the coherence of the 
propagation. In addition to these parameters, we shall also examine how source/receiver 
geometry, which affects propagation path, may influence the signal coherence. Since the 
propagation paths are related to the background sound speed profde, analysis will also be 
performed with other SSPs. 
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1. Turbulence Strength 

As noted in Chapter II, the spectral density of turbulence, given by Eq. (65), has 
an overall scaling factor, , which we refer to as the turbulent scaling parameter. Any 

desired RMS perturbation may be obtained by adjusting this scaling parameter. In this 
work, the turbulent scaling parameter has been fixed at five specific values, considered to 
range over realistic turbulence perturbation scales from very weak to very strong. One 
value was considered a reference value, two values were examined below reference, Ref- 
1 and Ref-2, and two values were examined above reference, Ref+1 and Ref+2. Sample 
values for the RMS perturbation for each of the turbulent strength cases is provided in 
Table 2. 


Case 

Ref-2 

Ref-1 

Reference 

Ref+1 

Ref+2 

RMS 

0.112 

0.224 

0.448 

0.897 

1.793 


Table 2. RMS Perturbation of Sound Speed (m/s) for Different Turbulent Strengths 

2. Turbulence Outer Length Scale 

The outer length cutoff scale, k ,, in Eq. (65) serves as a low wavenumber cutoff, 
limiting contributions from the turbulent perturbation to wavenumber values roughly 
greater than k,. Since 



(78) 


then \ is a wavelength threshold that limits the spatial scale of the turbulent fluctuations. 


The weak scattering theory dictates that the scattering of a propagating signal 
results mainly from perturbations on the scale of the Eresnel radius. The Eresnel radius, 
Rf , is given by (Duda et al, 1998) 

, (79) 
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where r is the transmission range and A, is the wavelength of the signal. For a signal with 
a frequency of 4 kHz, the wavelength is given by 


;l = 


/ 


1500m/sec 

4000Hz 


0.375 m. 


(80) 


Thus for a transmission range of 500 m, the Fresnel radius is 

^/, 5 oo«,,i 6 mz = V500m*0.375m □ 13.7m . (81) 

The corresponding “Fresnel wavenumber” is then about 0.46 m"'. Table 3 shows the 
Fresnel radius and Fresnel wavenumber values for frequencies examined in this research. 


Frequency (kHz) 

Wavelength (m) 

Fresnel Radius (m) 

Fresnel Wavenumber (1/m) 

4 

0.375 

13.693 

0.459 

8 

0.188 

9.682 

0.649 

16 

0.094 

6.847 

0.918 

20 

0.075 

6.124 

1.026 


Table 3. Fresnel Radius and Fresnel Wavenumber Values for Different Frequencies at a 
Range of 500 m 

Though the values for the Fresnel radius vary between 6 and 14 m, 10 m is used 
as the reference case for the turbulence outer length scale for all frequencies in the 
research. The other values examined in this work are 2 m and 5 m, as Ref-2 and Ref-l, 
respectively, and 25 m and 50 m, as Ref+1 and Ref+2, respectively. As the scale value 
was increased from low values, the most distinguishable change occurred around 10 m, 
which is in the vicinity of the Fresnel radius. So the Fresnel radius may be considered as 
a critical threshold for the values of outer scale regarding the disturbance of a signal by 
turbulence. In other words, since we assumed 10 m as the reference value for the 
turbulence outer length scale, the correlation curves should begin to exhibit a 
decorrelation at this value. As we increase this scale parameter to values greater than the 
Fresnel radius, then the correlation curves may change slightly. But as we decrease this 
parameter to values less than the Fresnel radius, a dramatic change in coherence should 
occur, and the correlation should increase considerably. 
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Figures 12 to 14 show one-dimensional slices of the turbulent perturbation energy 
density spectrums in the wavenumber domain for different outer length scales. It is 
important to note that around 2 m'\ the amplitudes of the spectra are very similar. 


Spectral Density of Turbulent Sound Speed Perturbation (1-D slice) 



Figure 12. One-dimensional energy density spectra for reference outer length scale (10 
m) and values below (2m and 5m). Vertical lines are the Fresnel radius values for 4, 8, 
16, and 20 kHz, from left to right, respectively. 
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Spectral Density of Turbulent Sound Speed Perturbation (1-D slice) 



Figure 13. One-dimensional energy density speetra for reference outer length scale (10 
m) and values above (25m and 50m). Vertical lines are the Fresnel radius values for 4, 8, 
16, and 20 kHz, from left to right, respectively. 


Spectral Density of Turbulent Sound Speed Perturbation (1 -D slice) 



Figure 14. One-dimensional energy density spectra for reference outer length scale (10 
m) and values above (25m and 50m) on a smaller scale. Vertical lines are the Fresnel 
radius values for 4, 8, 16, and 20 kHz, from left to right, respectively. 
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Another important point to make about the influence of the turbulent outer length 
scale is that changing this value changes the RMS sound speed perturbation value. For 
the reference turbulent strength scaling, sample RMS values for the five different outer 
length cut-off scales are shown in Table 4. 


Case 

Ref-2 

Ref-1 

Reference 

Ref+1 

Ref+2 

RMS 

0.217 

0.338 

0.448 

0.627 

0.799 


Table 4. RMS Perturbation of Sound Speed (m/s) for Different Cases of Turbulence 
Outer Length Scale Parameter 

3. Sound Speed Profiles 

During the ASIAEX, 54 CTD casts were made by the Chinese research vessel 
Shiyan-3 from 2 to 6 June, 2001. The reference SSP is taken as the average value of these 
measurements. In addition to the reference sound speed profile, three other individual 
sound speed profiles were examined in this analysis. One was chosen as being reasonably 
close to the mean (reference) profile, one was chosen as being the most deviated from 
mean with a deeper thermocline, and the third was chosen as being the most deviated 
from mean with the uppermost thermocline. These three profiles are referred to as 
Closest to Mean (CTM), Most Deviated Down (MDD), and Most Deviated Up (MDU), 
respectively. Figure 11 shows how the four sound speed profiles compare. 

Analysis of the influence of background sound speed profile is more complex 
than the other environmental parameters because the variability introduced in the signal 
coherence is not easily predictable. Indeed, we initially assumed that the background 
profile would not significantly affect such short range coherence results. The analysis 
presented in the following chapter, however, does indicate a significant dependence on 
the background profile and the propagation path geometry. 

Previous work by Duda and Bowlin (1994) examined the influence of sound 
speed relative curvature on signal stability. However, this approach was examined in the 
context of long-range, deep-ocean environments and assumed to hold for multi-cycle ray 
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propagation. Therefore, it is unclear if a similar measure of the background sound speed 
profiles in this environment will relate in any way to signal coherence along specific ray 
path geometries. 


Sound Speed Profiles 



Figure 11. Sound Speed Profiles Used in the Modeling Analysis 
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IV. POST PROCESSING AND ANALYSIS 


As explained in Chapter II, the vertieal spatial eorrelation is used to examine the 
signal degradation. If signal degradation is not present, then the eorrelation is going to be 
very elose to unity. A perfeet eorrelation value of 1 is the expeeted result of an 
autoeorrelation (the eorrelation of the signal with itself). All other eombinations for 
eorrelation are normalized aeeording to the autoeorrelation result, and thus they are less 
than or equal to 1. 

The eorrelations are ealeulated as a funetion of the “depth 
differenee/wavelength.” The eenter frequeney of the signal and Cg = 1500m/s are used 
for the wavelength estimation. 

For eaeh parameter, the referenee ease is shown first followed by the various 
results for other values of that parameter. Therefore, some figures are ineluded several 
times to allow for eonvenient eomparison. “Referenee Case” means that all parameters 
are adjusted to their referenee values. 

During the analysis of the turbulenee outer length seale, sound speed profde, and 
souree depth parameters, the referenee turbulenee strength didn’t ereate enough 
variability for 4 kHz and 8 kHz. Therefore, during the analysis of these parameters, the 
referenee seale for these frequeneies is ehanged. Speeifieally, the Ref+2 turbulenee 
strength parameter is used for the 4 kHz analysis, and the Ref+1 turbulenee strength 
parameter is used for 8 kHz analysis. All graphs exhibit four different eorrelation plots 
for the different reeeiver sub-arrays, with the uppermost plot eorresponding to the 
shallowest sub-array and the lowermost eorresponding to the deepest. 

A. TURBULENCE STRENGTH 

The turbulenee strength is a parameter that direetly affeets the scale of turbulence 
in the environment. As expected, the influence of this parameter is found to depend upon 
frequency of the propagation. Examining Figures 15 through 18, it can be easily seen that 
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as the frequency is increased, the same turbulence strength produces significantly 
more variability in the propagation. This effect is most noticeable when the frequency 
is doubled (from 4 kHz to 8 kHz, then 8 kHz to 16 kHz), but is not as evident above 16 
kHz in this analysis. Figures 17 and 18 show that 20 kHz propagation has almost the 
same vertical coherence as 16 kHz. 

Another important result to note is that the coherence for the deepest sub-array is 
a lot higher than the other sub-arrays for almost all cases of varying frequency or strength 
parameter. This may be related to the background sound speed profile with which the 
acoustic path interacts. This will be examined further in a subsequent section using a 
simple ray tracing model. 

Also note that these figures indicate little effect of the turbulence at 4 kHz when 
the reference strength parameter is employed. As previously noted, all subsequent 
analysis at 4 kHz will be performed using the highest (Ref+2) turbulence strength 
parameter. Similarly at 8 kHz, all subsequent analysis will be performed using the 
second highest (Ref+1) turbulence strength parameter. 
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(a) 



(b) 


(c) 



(d) (e) 

Figure 15. Effect of Turbulence Strength Parameter at 4 kHz. 

(a) Turbulence Strength = Reference, (b) Turbulence Strength = Ref-2, (c) Turbulence 
Strength = Ref-1, (d) Turbulence Strength = Ref+1, (e) Turbulence Strength = Ref+2 
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(a) 




(d) (e) 

Figure 16. Effect of Turbulence Strength Parameter at 8 kHz. 

(a) Turbulence Strength = Reference, (b) Turbulence Strength = Ref-2, (c) Turbulence 
Strength = Ref-1, (d) Turbulence Strength = Ref+1, (e) Turbulence Strength = Ref+2 
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(a) 



(b) 


(c) 



(d) (e) 

Figure 17. Effect of Turbulence Strength Parameter at 16 kHz. 

(a) Turbulence Strength = Reference, (b) Turbulence Strength = Ref-2, (c) Turbulence 
Strength = Ref-1, (d) Turbulence Strength = Ref+1, (e) Turbulence Strength = Ref+2 
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(d) (e) 

Figure 18. Effect of Turbulence Strength Parameter at 20 kHz. 

(a) Turbulence Strength = Reference, (b) Turbulence Strength = Ref-2, (c) Turbulence 
Strength = Ref-1, (d) Turbulence Strength = Ref+1, (e) Turbulence Strength = Ref+2 
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B. TURBULENCE OUTER LENGTH SCALE 


The analysis of the turbulence outer length scale confirmed that acoustic 
propagation through a turbulent environment is most sensitive to structures on the scale 
of the Fresnel radius. If the range is constant, the Fresnel radius is a function of frequency 
only. In this model, the range is kept constant and the frequency range of 4-20 kHz is 
used. As a result the Fresnel radius values vary from 13.7 m down to 6.1 m. Rather than 
changing the cutoff scale for each frequency, kt = 0.628 1/m (outer length scale of /l(= 10 
m) is used as a reference value, consistent with previous work (Duda et ah, 1998). Note 
that this falls within the range of Fresnel scales defined. Thus, the influence of the outer 
length scale is considered here. 

Figures 19 to 22 show the effect of changes in cutoff scale at different 
frequencies. The vertical correlation is found to exhibit the greatest degradation (most 
significant scattering) for the largest outer length scale value of 50 m. This result is 
verified for all frequencies. 

A turbulence outer length scale of 5 m barely produced any visible signal 
degradation, and a cutoff scale of 2 m produced no observable signal degradation for all 
frequencies. For these cutoff scale values, the turbulence spectrum does not generate 
enough perturbation at the Fresnel scale for the frequencies considered. However, as the 
cutoff scale is increased to 10 m, the part of the spectrum generating perturbations at the 
Fresnel scale starts to increase, and the signal correlation tends to decrease. This result is 
then consistent with the expectation of the dominant effect of scattering at the Fresnel 
scale. 

For 4 kHz, the Fresnel scale is 13.7 m. When the reference cutoff scale of 10 m is 
applied, the perturbations are below the Fresnel scale. Therefore, a cutoff scale value of 
25 m has a remarkable effect for this frequency. The same transition may also be 
observed for 8 kHz, since the Fresnel scale value for this frequency is 9.7 m, which is 
slightly below the reference cutoff scale. The Fresnel scale values for 16 and 20 kHz are 
6.8 and 6.1 m, respectively. When the reference cutoff scale is applied, the perturbations 
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are already beyond the eritieal threshold. Therefore, inereasing the eutoff seale to 25 m 
for these frequeneies does not eause a remarkable degradation in signal eoherenee. 

Comparison of the analysis of 8 kHz versus 16 and 20 kHz verifies that the 
dissipation of this effeet does not oeeur right beyond the threshold, but several meters 
beyond the threshold. When a 10 m eut-off seale was applied for 8 kHz, the perturbations 
are only 0.3 m beyond the threshold. However, these same perturbations are 3.2 and 3.9 
m beyond the threshold for 16 and 20 kHz, respeetively. The impaet of this ehange of 
seale on the level of perturbation at the Fresnel seale is readily observed in the speetral 
plots, displayed previously in Figures 13 and 14. Though an exaet value was not 
ealeulated, we eoneluded that if the perturbations are several meters greater than the 
Fresnel seale, further inerease of the eut-off seale doesn’t signifieantly affeet the signal 
eoherenee. Figures 19 to 22 elearly display the verifieation of these results. 

Although larger eutoff seales eould have been used for further analysis, the 
referenee value of 10 m was ehosen for the remainder of this thesis. The deeision was 
made so that the eutoff seale was eonsistent with previous researeh, and beeause this 
seale did produee notieeable signal degradation. Therefore, the infiuenee of other 
parameters eould be observed at this eutoff seale. 
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(a) 



(b) (c) 



(d) (e) 

Figure 19. Effect of Turbulence Outer Length Scale at 4 kHz. 

(a) Turbulence Scale = 10 m (Reference Case), (b) Turbulence Scale = 2 m, (c) 
Turbulence Scale = 5 m, (d) Turbulence Scale = 25 m, (e) Turbulence Scale = 50 m 
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(a) 



(b) 


(c) 



(d) (e) 

Figure 20. Effect of Turbulence Outer Length Scale at 8 kHz. 

(a) Turbulence Scale = 10 m (Reference Case), (b) Turbulence Scale = 2 m, (c) 
Turbulence Scale = 5 m, (d) Turbulence Scale = 25 m, (e) Turbulence Scale = 50 m 
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(b) 


(c) 



(d) (e) 

Figure 21. Effect of Turbulence Outer Length Scale at 16 kHz. 

(a) Turbulence Scale = 10 m (Reference Case), (b) Turbulence Scale = 2 m, (c) 
Turbulence Scale = 5 m, (d) Turbulence Scale = 25 m, (e) Turbulence Scale = 50 m 
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(a) 



(b) 


(c) 



(d) (e) 

Figure 22. Effect of Turbulence Outer Length Scale at 20 kHz. 

(a) Turbulence Scale = 10 m (Reference Case), (b) Turbulence Scale = 2 m, (c) 
Turbulence Scale = 5 m, (d) Turbulence Scale = 25 m, (e) Turbulence Scale = 50 m 
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c. 


CURVATURE ANALYSIS 


Previous work by Duda and Bowlin (1994) examined the stability of long-range 
aeoustie ray propagation through internal wave fluetuations in deep water. They found 
that the stability of the rays propagating through various depths of the waveguide were 
related to mean values of loglf/l, where 


U = c 






(93) 


and c = c(z) is the depth-dependent sound speed profde. More reeent work by Brown, et 
al. (2004, in press) also showed theoretieally how long-range ray stability may be 
dominated by features of the baekground sound speed profde. 

While both of the aforementioned analyses were eoneemed with long-range, deep 
water propagation, the analysis presented here suggests a potentially similar influenee of 
the baekground sound speed profile. Therefore, we shall now examine a similar measure 
of the direet path propagation in the shallow water environment presented here. 
Speeifieally, the depth-dependent parameter U , defined above, will be eomputed for the 
four profiles examined. A simple ray-traee model will then be employed to eompute ray 
trajeetories through the unperturbed (no turbulenee) baekground profiles and the 
subsequent average of /og|f/| along the ray path will be determined. While Duda and 
Bowlin (1994) ineluded their internal wave perturbations in the ealeulations of U, the 
turbulenee fluetuations are not expeeted to signifieantly affeet the average of /og|f/| 
along the direet-path ray trajeetories. 

In order to eompute U, the diseretely sampled sound speed profiles provided by 
the measured data sets must be fitted by some twiee differentiable funetion. Cubie 
splines are employed here. Although Duda and Bowlin (1994) refer to possible spurious 
values from sueh eubie polynomials, this tends to be a more signifieant problem when the 
spaeing between measured sound speed values is not uniform in depth. In the ease of the 
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ASIAEX data sets used here, the depth sampling is uniform, so the use of cubic 
polynomials is expected to be more stable. 

As an example of the data fit achieved by the cubic spline technique. Figure 23 
displays the CTM profile data (as discrete points) compared to data interpolated from the 
cubic polynomial. This comparison shows the density of measured data points relative to 
the observed curvature in the profile, and justifies the use of the cubic splines to fit the 
data. A direct calculation of log|l7| is then possible. Figure 24 provides the results of 
these calculations for all four sound speed profiles considered. 



Figure 23. Measured Sound speed Data Compared to Cubic Spline Fit for CTM Profile. 
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Log |U| (MOD) 


Log |U| (MDU) 


Figure 24. Profiles of log|t/| for the four environmental profiles examined: reference 
(upper left); CTM (upper right); MDD (lower left); MDU (lower right). 


To compute ray trajectories, a simple 4*-order Runge-Kutta algorithm was 
employed. Such an algorithm requires a continuous sampling of both the sound speed 
profile and its first derivative, both of which are provided by the cubic spline 
interpolation. Although the code can accurately handle smooth boundary reflections, such 
ray paths were not of interest here, so only those rays without boundary reflections were 
saved for analysis. 

At each range step along each ray path, the value of log|[/| was then evaluated. 
By simply summing all such values and dividing by the number of range steps in the 
calculation out to 500m, an along-ray average of loglt/l was computed. In Figure (4), 
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sample plots of ray paths in the reference environment for three different source depths 
(30m, 50m, and 70m) are displayed. It is important to note that, in order to compute a 
well-sampled profile of the along-ray mean of log|f/|, a dense sampling of ray arrival 

depths at the final range (500m) is needed. Because the structure of the background 
sound speed profiles tended to cause a divergence of ray paths near the center of the 
water column (refer to 50m source depth ray paths in Figure 25), various techniques were 
used to increase the density of ray launch angles near 0“, and the minimum and maximum 
values of launch angles were varied to insure arrival depths between the depths of interest 
(about 20m-80m). Results of this analysis are also displayed in Figure 25 for the same 
environment and source depths, but using an altered code that increases the density of ray 
paths for adequate sampling of the final profile of the along-ray mean of log|f/|. 
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Figure 25. Sample ray paths for the reference environment and the corresponding profiles 
of the along-ray mean of log|t/| for three source depths: 30m (upper plots); 50m (middle 
plots); 70m (lower plots). 

The profiles of the along-ray mean of log|l7| were computed for all background 

sound speed profde environments and source depths considered. Once these profiles 
were determined, a final average value over the depth of the sub-arrays was then 
computed. These unique values are tabulated and presented in the following sections to 
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compare with the observed dependence of signal correlation on propagation paths 
through the various environments. 

D. SOURCE DEPTH 

The depth of the souree will have a direet influenee on the propagation paths 
between the source and reeeiver elements as dietated by the sound speed profde. The 
reference SSP is used for the analysis of source depth influence. The reference SSP is 
approximately isospeed over the depth intervals 0-20 m and 80-100 m, and as a 
eonsequence doesn’t exhibit a remarkable variability at these depths. There is a 
thermocline between the depths 20-80 m, and so propagation through this region 
undergoes significant refraction. 

As the source depth was increased from 30 to 40 and then 50 m, an observable 
pattern existed in the change of signal eoherence for all frequeneies. When the souree 
was placed at 30 m, the signal mostly propagated in the upper, isospeed part of the 
waveguide. These paths appear coincident with little signal degradation. As the depth of 
the souree is lowered to 40 m, the signal partly propagates through the thermocline. This 
eaused a slight deerease in the correlation for all depths and frequeneies. When the 
source depth was ehanged to 50 m, the sub-arrays at 45 and 55 m reeeived signals whieh 
traversed along the thermoeline. The signals recorded at these sub-arrays exhibit a 
remarkable drop in the signal eoherenee. Verification of these results is shown in Figures 
26 to 29 (parts a, b and e of figures). 

When the source was lowered to 60 m, signals at the upper sub-arrays, 35 and 45 
m, observed significantly more variability compared to the lower sub-arrays, consistent 
with higher variability along those paths propagating through regions of higher sound 
speed gradient. Plaeing the source at a depth of 70 m didn’t eause a noteworthy 
degradation of correlation at any receiver depths for all frequencies tested in this 
researeh. 

Due to this apparent eorrelation between background sound speed gradient, 
propagation paths, and signal eoherenee, the analysis of the variable source depth effeet 
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includes the effeet of SSP eurvature, as described in the previous seetion. The effect of 
SSP curvature was not consistent for all results of souree depths. However, considering 
that this effect was originally derived for deep and long-range propagation, the 
exeeptions may be due to the shallow water, short-range nature of the problem. 

The sub-array averages of eorrelation values are calculated to provide a 
eomparison between the signal coherenee and SSP curvature. All of the eorrelation 
eurves generally exhibited a deerease in eoherenee over the first few wavelengths, and a 
slightly inereasing coherence pattern near the last wavelength. The stabilized region of 
the correlation curve is extracted after removing the correlation values corresponding to 
the first five and the last wavelengths. These stabilized correlation eurve values were then 
used to generate an average eoherenee value for each sub-array. 

When there was a noticeable relative differenee between the eorrelation values of 
receivers, the relative ehange of the SSP eurvature parameter was eonsistent. For the 
souree depths of 30 and 40 m, a significant decorrelation was not noted at any depths and 
frequeneies. Therefore, the infiuenee of SSP curvature isn’t analyzed for these souree 
depths. The sub-array averages of along-ray means of log|f/| for the reference sound 
speed profile and depths of 50, 60, and 70 m are exhibited in Table 5. The sub-array 
averages of eorrelation values for souree depths of 50, 60, 70 m are exhibited in Tables 6 
to 8. 


Receiver 

50 m 

60 m 

70 m 

35-40 m 

3.54 

3.33 

3.26 

45-50 m 

3.43 

3.13 

2.97 

55-60 m 

3.06 

2.82 

3.00 

65-70 m 

2.76 

3.00 

3.03 


Table 5. Sub-Array Averages of Along-Ray Means of log|f/|. 
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Receiver''''*^^ 

4 kHz 

8 kHz 

16 kHz 

20 kHz 

35-40 m 

0.91 

0.88 

0.94 

0.92 

45-50 m 

0.75 

0.82 

0.85 

0.89 

55-60 m 

0.85 

0.88 

0.89 

0.88 

65-70 m 

0.95 

0.97 

0.98 

0.98 


Table 6. Sub-Array Averages of Correlation Values for Souree Depth = 50 meters. 


Receiver''''''''''''^^ 

4 kHz 

8 kHz 

16 kHz 

20 kHz 

35-40 m 

0.86 

0.86 

0.90 

0.85 

45-50 m 

0.78 

0.76 

0.84 

0.89 

55-60 m 

0.94 

0.94 

0.97 

0.93 

65-70 m 

0.94 

0.95 

0.97 

0.97 


Table 7. Sub-Array Averages of Correlation Values for Souree Depth = 60 meters. 


Receiver''''**'''*^^ 

4 kHz 

8 kHz 

16 kHz 

20 kHz 

35-40 m 

0.84 

0.89 

0.92 

0.87 

45-50 m 

0.87 

0.89 

0.89 

0.93 

55-60 m 

0.92 

0.86 

0.93 

0.89 

65-70 m 

0.95 

0.93 

0.93 

0.93 


Table 8. Sub-Array Averages of Correlation Values for Source Depth = 70 meters. 


As the average along-ray mean values of log|f/| corresponding to each sub-array 
increase, the signal coherence is expected to decrease. This negative effect of SSP 
curvature may be observed by cross-examining Tables 5 to 8. There are exceptions of this 
effect, which again may be due in part to the shallow water, short-range nature of the 
problem. 
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(a) 



(b) (c) 



(d) (e) 


Figure 26. Effect of Source Depth Variation at 4 kHz. 

(a) Source Depth = 50 m (Reference Case), (b) Source Depth = 30 m, (c) Source Depth = 
40 m, (d) Source Depth = 60 m, (e) Source Depth = 70 m 
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(d) (e) 


Figure 27. Effect of Source Depth Variation at 8 kHz. 

(a) Source Depth = 50 m (Reference Case), (b) Source Depth = 30 m, (c) Source Depth = 
40 m, (d) Source Depth = 60 m, (e) Source Depth = 70 m 
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(b) 


(c) 



(d) (e) 

Figure 28. Effect of Source Depth Variation at 16 kHz. 

(a) Source Depth = 50 m (Reference Case), (b) Source Depth = 30 m, (c) Source Depth = 
40 m, (d) Source Depth = 60 m, (e) Source Depth = 70 m 
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(a) 



(b) 


(c) 



(d) (e) 

Figure 29. Effect of Source Depth Variation at 20 kHz. 

(a) Source Depth = 50 m (Reference Case), (b) Source Depth = 30 m, (c) Source Depth = 
40 m, (d) Source Depth = 60 m, (e) Source Depth = 70 m 
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E. SOUND-SPEED PROFILE 


Four different baekground sound speed profiles were examined in this researeh. 
Details of these profiles were deseribed in Chapter III. During the following analysis of 
the effeet of baekground SSP, the souree depth was fixed at 50 m. 

As noted previously, the referenee SSP is nearly isospeed in the intervals 0-20 m 
and 80-100 m, and eontains a thermoeline between the depths 20-80 m. The signal 
eoherenee was already observed to suffer the most degradation when the signal 
propagates through the thermoeline. Sub-arrays loeated at 45-50 and 55-60 m were found 
to reeeive the most degraded signals when the referenee SSP is used. This observation is 
verified for all frequeneies and the results are exhibited in Figures 30 to 33 (upper left 
parts of all figures). 

The CTM sound speed profile has depth properties similar to the referenee SSP. 
However, the CTM profile has more variability in the thermoeline region. This inerease 
of variability in sound speed is found to produee even more degradation in signal 
eoherenee. A deerease in eorrelation at all depths is observed. It must be noted that the 
signals on the sub-arrays loeated at 45-50 and 55-60 m still exhibit the greatest 
degradation eompared to other sub-arrays, presumably sinee they are loeated in the 
middle of thermoeline region. These results may be eonfirmed by looking at Figures 30 
to 33 (upper right parts of all figures). 

The MDD sound speed profile is mostly isospeed over the depths 0-65 m (with 
the exeeption of a single “jump” around 30 m). Therefore, based on the previous 
findings, we may expeet that signals from a 50 m souree will not be as severely degraded 
as they are with the referenee profile, sinee most of the propagation oeeurs in this upper, 
isospeed portion of the profile. The analysis of the results using the MDD sound speed 
profile eonsistently verified this assumption and the results are displayed in Figures 30 to 
33 (lower left parts of all figures). 

The MDU sound speed profile has effeetively no upper isospeed region. The 
thermoeline exists from around 5 m down to around 65 m, below whieh the profile is 
nearly isospeed. The overall variability of the MDU sound speed profile appears less than 
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the CTM sound speed profile. Therefore, it may be antieipated that the results of these 
two sound speed profdes from a 50 m souree would be similar, with the MDU profde 
having a weaker effeet on signal eoherenee. The model verified these assumptions and 
the results are shown in Figures 30 to 33 (upper right parts of all figures). 

Similar to analysis of souree depth, effeet of SSP eurvature is also analyzed for 
ehanges in the baekground sound speed profile. The most eonsistent results were noted 
for the referenee SSP and the CTM sound speed profile. The sub-array averages of along- 
ray means of log|f/| for the 50 m souree depth in the various sound speed profiles are 
tabulated in Table 9. The sub-array averages of eorrelation values for the various sound 
speed profiles are displayed in Tables 10-11. 


Receiver''''*''''^^ 

Reference 

CTM 

MDD 

MDU 

35-40 m 

3.54 

4.60 

4.67 

4.62 

45-50 m 

3.43 

4.86 

5.30 

4.27 

55-60 m 

3.06 

4.05 

4.38 

4.26 

65-70 m 

2.76 

3.69 

3.87 

4.26 


Table 9. Sub-array averages of along-ray means of log|f/| for the source at 50 m depth. 


Receiver''''''''*'^^ 

4 kHz 

8 kHz 

16 kHz 

20 kHz 

35-40 m 

0.91 

0.88 

0.94 

0.91 

45-50 m 

0.75 

0.82 

0.85 

0.75 

55-60 m 

0.85 

0.88 

0.89 

0.85 

65-70 m 

0.95 

0.97 

0.98 

0.95 


Table 10. Sub-Array Averages of Correlation Values for the reference SSP. 
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Reeeiver''''''''*''^^ 

4 kHz 

8 kHz 

16 kHz 

20 kHz 

35-40 m 

0.94 

0.94 

0.94 

0.88 

45-50 m 

0.72 

0.71 

0.73 

0.78 

55-60 m 

0.67 

0.65 

0.74 

0.74 

65-70 m 

0.79 

0.85 

0.88 

0.90 


Table 11. Sub-Array Averages of Correlation Values for the CTM sound speed profile. 


As the average along-ray mean values of log|f/| eorresponding to eaeh sub-array 
inerease, the signal eoherenee may be expeeted to deerease. The results of the analysis 
presented here verify this expeetation with some exeeptions. The results may be best 
understood by eross-examining Tables 9 to 11. The exeeptions are possibly due to the 
short-range nature of the propagation examined in this work. 
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(a) 


(b) 



(c) (d) 

Figure 30. Effect of Background Sound-Speed Profile at 4 kHz. 

(a) SSP = Reference Case, (b) SSP = CTM, (c) SSP = MDD, (d) SSP = MDU 
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Figure 31. Effect of Background Sound-Speed Profde at 8 kHz. 

(a) SSP = Reference Case, (b) SSP = CTM, (c) SSP = MDD, (d) SSP = MDU 
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Figure 32. Effect of Background Sound-Speed Profde at 16 kHz. 

(a) SSP = Reference Case, (b) SSP = CTM, (c) SSP = MDD, (d) SSP = MDU 
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Figure 33. Effect of Background Sound-Speed Profde at 20 kHz. 

(a) SSP = Reference Case, (b) SSP = CTM, (c) SSP = MDD, (d) SSP = MDU 
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V. SUMMARY 


A. CONCLUSIONS 

The motivation of this research was to develop a better understanding of the 
influence of turbulence and environmental variability on broadband acoustic coherence. 
The MMPE model was used to simulate the effects of a turbulent environment on such 
broadband acoustic propagation. Realistic turbulence perturbations based on well 
established statistical models were incorporated into the MMPE model, making 
investigations of such effects on signal coherence easily accessible for a variety of 
environmental factors and source/receiver geometries. The background sound speed 
profiles and the model geometry were similar to a component of one of the East China 
Sea experiments conducted during ASIAEX. Effects of turbulence strength, turbulence 
cut-off scale, background sound speed profde, frequency and source/receiver geometries 
on short range, direct-path propagation were examined. 

Vertical correlation was used to quantify the coherence of the received signal at 
four vertical sub-arrays located at a fixed range from the source. Each case was simulated 
10 times to simple statistics (mean and standard deviation) of the signal decorrelation due 
to turbulence. 

It was observed that the turbulence strength, a parameter that directly affects 
the RMS magnitude of the turbulence perturbation, exhibited a positive correlation with 
signal degradation for all frequencies up to 20 kHz, as expected. Essentially, this simply 
states that the stronger the sound speed perturbations are, the stronger the acoustic signal 
degradation will be. 

The analysis of the turbulence outer length scale also confirmed the expected 
result that acoustic propagation through a turbulent environment is most sensitive to 
structures on the scale of the Eresnel radius. When the perturbations are already greater 
than the Eresnel scale, further increase in cut-off scale (even with increases in total RMS 
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perturbation strength) didn’t affect the signal coherence. An exact value was not 
calculated for the critical excess value beyond which the coherence is not affected 
remarkably. However, it was observed that if the turbulence cut-off scale is several 
meters greater than the Fresnel scale, then further increases of the cut-off scale do not 
produce a noteworthy effect on the signal coherence. 

Source-receiver depth and background SSP were found, somewhat surprisingly, 
to have a significant influence on the signal coherence. This was not expected because 
there was no direct relationship between the background profile and the structure of the 
turbulence, and because the propagation distances were so short that such influences 
were not expected to have enough time to “build up” along the path. It was confirmed 
that propagation through a thermocline has greater (negative) effect on signal coherence 
than the isospeed regions of the SSP. Thus, the signal stability appears to be strongly 
dependent upon the propagation path through the background profde (without 
turbulence). 

Due to these observations, an attempt was made to compare these results with 
similar findings of signal stability dependence in deep ocean, long-range propagation on 
sound speed profde curvature. The analysis of the SSP curvature performed here was 
partly consistent with the findings of the signal coherence. Most notably, when there 
were significant relative differences between the correlation values of the sub-arrays, the 
relative differences of the mean curvature parameter along the ray paths were consistent. 

B. RECOMMENDATIONS AND FUTURE WORK 

The findings of the significance of the structure of the background sound speed 
profile were rather surprising, but of considerable significance. Such findings could result 
in changes in the operation of various underwater acoustic systems. For example, tactical 
systems (e.g., mine-hunting UUVs) may be able to improve their detection capabilities by 
determining their optimal operating depths based on in situ measurements of the 
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background profile. Depth variable aeoustie eommunieation systems may also be able to 
optimize their performanee in a similar fashion. 

While the analysis of SSP eurvature performed in this thesis showed some 
eonsisteney with the signal eoherenee results, it is important to reeall that this measure 
was originally derived for deep and long-range propagation. It is possible, therefore, that 
another measure of the interaetion of the propagation paths with the baekground profile 
may be better suited for predietion of signal degradation. Future work may then examine 
other such measures. 

It should also be noted that this analysis assumed a uniform turbulenee 
perturbation field. The effeets of depth or range variability of the turbulenee is unknown, 
as is any obvious eorrelation with depth or range variability of the baekground profile. 
More sophistieated oeeanographie models of turbulenee would be needed to generate the 
neeessary realizations of turbulent fiuetuations in order to study these effeets. 

Finally, although this modeling effort was based on the experimental geometry 
employed during the ASIAEX studies in the East China Sea, the measured data sets did 
not inelude any direet measurements of turbulenee structures. It would be of great 
benefit in a future experiment to simultaneously measure both signal eoherenee and water 
eolumn fine-seale strueture. The validity of this numerieal approaeh may then be 
eonfirmed, and further numerieal studies eould then be performed with eonfidenee. Our 
understanding of the impaet of small-seale oeean fiuetuations on various aeoustie 
systems’ performanee eould then be enhaneed, and future system designs eould be 
improved. 
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